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1.0  INTRODUCTION 

In  the  past  decade,  experimental  and  theoretical  investigations 
at  a large  variety  of  industrial,  governmental,  and  university  research 
centers  have  revealed  an  unexpected  complexity  in  the  physics  of  hyper- 
sonic phenomena.  The  quest  for  new  simple  fundamentally  consistent  and 
universal  scaling  laws,  even  in  nonreacting  flows,  has  been  repeatedly 
frustrated.  Furthermore,  the  increasingly  mysterious  coupling  of 
turbulence  and  molecular  interactions  at  high  Reynolds  numbers  and 
high  Mach  numbers  has  introduced  new  dimensions  to  the  uncertainty 
ordinarily  associated  with  the  predictability  of  high-speed  flow. 

By  contrast  (as  a result  of  the  same  investigations),  the  capability 
of  the  systems  designer  for  detailing  the  requirements  of  specific  systems 
has  increased.  This,  in  turn,  enhances  the  likely  impact  on  the  current 
understanding  of  high-speed  phenomena  of  a carefully  planned  and  focused 
attempt  at  simulation.  Thus,  both  to  serve  the  interests  of  current 
system  problem-solvers  and  to  attack  anew  the  uncertain  physics  as  well, 
the  elements  of  model  similitude  deserve  a re-examination. 

With  this  motivation,  the  problem  of  slender  cone  simulation 
will  be  reviewed.  There  have  been  more  than  100  years  of  attention 
given  to  the  problem  of  gas  dynamic  similitude  (Ref.  1).  In  this  report, 
the  approaches  developed  by  Zierep  (Ref.  2)  and  Thompson  (Ref.  3)  are 
followed,  and  the  theoretical  analyses  of  Eschenroeder  (Ref.  A)  have 
been  used  extensively.  After  a brief  summary  of  theoretical  similitude, 
the  various  flow  regimes  will  be  cataloged.  Finally,  the  overall 
approach  in  this  report  will  be  applied  to  a specific  example  of  current 
interest. 
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2.0  SIMILITUDE:  A REVIEW 

The  variety  of  ways  in  which  a physical  law  can  be  expressed  is 
limited  by  the  requirement,  among  others,  that  its  form  be  independent 
of  a change  in  dimensions.  Thus,  even  when  the  law  is  not  known,  cer- 
tain restrictions  can  be  defined.  Consider,  therefore,  that  a list  of 
important  dimensional  parameters  for  a given  observable  (X^)  can  be 
defined.  When  this  list  is  ordered,  one  has  (including  X^): 


The  law  characterizing  these  parameters  can  always  be  put  in  the  form 

G (X.,  X_,  ...  X,,  ...  X ) = 0 
i z x n 

Usually,  there  are  only  p of  these  parameters  (p  <_  n),  which  are  dimen- 
sionally independent;  that  is,  only  (n  - p)  parameters  have  dimensions 
which  can  be  found  as  products  of  the  dimensions  of  one  or  more  of  the 
remaining  p parameters.  Further,  if  the  number  of  fundamental  dimensions 
(such  as  mass,  length,  time,  charge,  and  temperature)  relevant  to  the 
problem  is  q,  then  p must  be  equal  to  or  less  than  q.  Let's  rearrange 
our  list  of  X’s  into  two  groups  so  that  the  dimensionally  independent 
parameters  are  listed  first: 

■ 

x. , X.,  ...  x , : x x 

12  p . p+1  , n 

The  dimensions  of  the  (n  - p)  parameters  in  the  second  group  (j  p + 1 ) , 

Pd  ■ N 6'  Pd  '2  -•  Pd  6p 

can  now  be  expressed  in  terms  of  the  first  group  and  a dimensionless 
parameter  in  each  case  results  when  the  following  prescription  is  used: 


X. 
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From  this,  an  alternative  and  equivalent  form  of  the  physical  law  is 

b <V  V •••  Vp>  ’ 0 

as  the  law  now  allows  only  (n  - p)  independent  variables. 

If 'the  form  of  the  law  is  assumed  to  be  known,  then  there  are  at 
least  two  resulting  implications.  First,  it  is  implied  that  a decision 
has  already  been  made  concerning  the  important  physical  parameters. 
Secondly,  an  adequate  set  of  boundary  conditions  on  any  relevant  differen- 
tial equations  must  exist  and  be  expressable  in  terms  of  dimensionless 
dependent  and  independent  variables  and  using  nondimens ional  constant 
coefficients. 

Thus,  the  similarity  analysis  allows  a reduction  in  and  an  inter- 
dependence among  the  significant  physical  parameters.  In  addition,  for 
gas  dynamic  problems,  the  similarity  solution  is  a prerequisite  for  the 
development  of  model  rules  since  it  is  desired  to  transform  the  data  of 
flow  determined  from  one  body  to  another  body.  Sometimes,  it  is  even 
possible  through  the  modeling  transformation  to  obtain  insight  not  other- 
wise available  by  treating  an  ostensibly  unrelated  problem  whose  equations 
can  be  placed  in  a form  similar  to  the  gas  dynamic  equations. 

3.0  THE  SHARP  COIME  IN  AN  IDEAL  GAS 

Now,  one  can  apply*  these  considerations  to  the  general  problem  of 

high-speed  gas  flow  around  sharp  cones.  That  is,  for  cone  angle  (8), 

2 

it  is  assumed  that  sin  0 <<  1.  After  first  discussing  flow  without 
chemistry,  the  treatment  will  be  generalized  to  include  reacting  species. 

3.1  EQUILIBRIUM  FLOW 

Let's  assume  that  the  overall  field  is  describable  (for  an  ideal 
gas)  by  eighteen  parameters  (n  = 18),  allowing  for  all  fluid  dynamical 
possibilities,  as  shown  below: 
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1 . 

X*  ^ 

’ l 

I 

11 . o; 

surface  tension 

2. 

y! 

> spatial  coordinates 

12.  A; 

conductivity 

3. 

z-  ' 
z*  > 

I 

13.  c ; 

specific  heat. 

4. 

w; 

characteristic  velocity 

r 

constant  pressure 

5. 

p; 

pressure 

14.  C ; 

V 

specific  heat, 
constant  volume 

6. 

pi 

density 

15.  t; 

time 

7. 

characteristic  length 

16.  Q; 

energy  source  or  sink 

8. 

g; 

acceleration  due  to 
gravity 

17.  v; 

viscosity 

9. 

e; 

coefficient  of  expansion 

18.  a; 

heat  transfer  rate 

10. 

AT; 

temperature  difference 

Take  the  value  p = 4 (as  this  is  the  maximum  number  corresponding  to 

mass  length,  time,  and  temperature  as 

basic  units)  and  use 

X1 

= P» 

X_  = C , X0  = v,  and  X.  = l. 
1 p 3 4 

This  leaves  (n  - p = 14)  fourteen 

independent  variables: 

1. 

x;  x 

:/& 

9. 

2 

a;  Jlpw  /o  or  We 

2. 

y;  y/£ 

10. 

B;  gi^SAT/v^  or  Gr 

3. 

Cj 

C /C  or  y 

n 

V 

P v 

11. 

t;  ! C p/At  or  Fo 

4. 

w;  wi/v  or  Re 

P 

5. 

z;  z 

\/l 

12. 

a;  a/wpCp  or  St 

6. 

A;  vpC  /A  or  Pr 
P 

13. 

AT;  w2/C  AT  or  Ec 

P 

7. 

2 2 
p;  w p/yp  or  M 

14. 

Q;  Q/CpAT  or  D 

8. 

2 

g;  w /g£  or  Fr 

where 

Re 

= Reynolds  number 

Gr  = Grashof  number 

Pr 

= Prandtl  number 

Fo  = Fourier  number 

M 

= Mach  number 

St  = Stanton  number 

Fr 

= Froude  number 

Ec  = Eckert  number 

We 

= Weber  number 

D = Damkohler  number 
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Thus,  a general  solution  for  this  problem  can  be  expressed  in  terms  of 
the  well-known  similarity  parameters: 

b (r/l;  y.  Re,  Pr,  M,  Fr,  We,  Gr,  Fo,  St,  Ec,  D)  = 0 

When  this  result  is  applied  to  the  general  problem  of  modeling 
hypersonic  flow,  it  is  clear  that  complete  simulation  of  the  flow  is 
extremely  complicated.  Lrt’s  look  at  three  limiting  cases.  For 
frictionless  flow  and  no  important  unsteady,  gravitational  or  ther- 
modynamic effects. 


b^r/Jl;  y,  M)  = 0 

and  simulation  is  trivial.  If  viscous  effects  are  important  as  well, 

b2(r/£;  y,  M,  Re)  = 0 

This  implies  that,  ignoring  thermodynamics,  the  modeling  can  be  accom- 
plished for  viscous  flows  by 

b2a(r/A,  Mn,  Rem)  = 0 


Let 


_ n i m 

Z.  = M /Re 

where  n and  m are  determined  by  the  parameter  of  interest.  For  example, 

n = m = 1 corresponds  to  Z = Kundsen  number;  here,  Z is  the  well-known 

characterization  of  the  transition  from  continuum  to  free  molecular  flow. 

For  n = 3 and  M = 1/2,  Z characterizes  the  influence  of  shock  wave-boundary 

3 

5^  M 

layer  interference  (since  “v — ^ -,_).  Figure  1 shows  these  possible 

1 Re  ' L 

regimes.  Adding  unsteadiness,  the  equation  for  b(Y^)  becomes 
b3(r/t,  Y,  M,  Re,  Fo)  = 0 
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In  this  last  case,  the  modeling  and  approximating  problem  is  obviously 
more  complicated. 

100.0 

60.0 


10.0 

6.0 

M 


1.0 
0.6 

10"1  101  It P 10?  107  109 

Re 

Figure  1.  Boundaries  of  flow  regimes  in  the  M-Re  plane  {from  Ref.  2). 

3.2  NON-EQUILIBRIUM  EFFECTS 

When  molecular  and/or  atomic  interactions  require  a period  of  time 
in  order  to  reach  equilibrium  (t^)  which  is  comparable  to  local  flow 
field  time  (tf  O^/w)),  then  it  becomes  necessary  to  introduce  a second 
time  variable.  This  results  -in  a new  time  similarity  parameter: 
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Either  t,  >>  t or  tc  « t will  make  D irrelevant.  However,  if  D is 
f r f r 

important,  then,  in  order  to  scale  both  D and  Re  in  modeling,  since 


D 


and  Re  = 


w£ 

v 


the  different  dependence  on  w requires  that  one  be  able  to  vary  t and/or 
v.  Thus,  an  important  limitation  is  imposed  on  these  parameters. 


If  there  are  several  processes  with  relaxation  times  (t^),  then 
each  one  has  a parameter  (D^)  given  by 


Di  = r 


r . 
1 


When  the  slender  cone  in  air  is  considered  (Ref.  4),  for  example,  the 

til  til 

i two-body  forward  reaction  for  the  j species  has  the  form 

(2)  _ P~*Vii  kfi(2)(Io) 
fij  w)li 


where  is  the  mean  molecular  weight  of  the  cold,  undisturbed,  flight 
medium,  Tq  is  a typical  flow  field  temperature,  . are  the  stoichiometric 
coefficients,  are  the  rate  constants,  and  the  other  symbols  have 
their  usual  meaning.  Similarly,  for  three-body  reactions: 


fij 


(3)  _ P°°  lvHkfj(  )(V 


w3l> 


In  the  above,  the  gas  dynamic  equations  have  been  used,  the  diffusion 

coefficient  is  assumed  to  be  proportional  to  Ta/p,  and  viscosity  is 

proportional  to  (T)p  where  a and  & are  unknown  powers.  From  this,  the 

simultaneous  preservation  of  both  D's  would  not  generally  be  possible 

2 

since  the  requirements  of  p^S,  and  £ are  contradictory. 
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Obviously,  the  long  list  of  similarity  parameters  must  be  reduced 
if  modeling  is  to  be  possible.  This  might  be  accomplished  by  a space- 
time  transformation  so  that  the  general  three-dimensional  steady  flow 
problem  becomes  a problem  in  one-dimensional  unsteady  flow.  This  might 
be  accomplished,  alternatively,  by  determining  and  exploiting  the 
flexibility  of  neglecting  certain  effects  in  the  various  flight  regimes. 
Furthermore,  the  similarity  problem  might  be  simplified  by  exactly 
reproducing  certain  aspects  of  the  prototype  problem  so  that  the  range 
of  modeling  to  be  attempted  is  reduced. 

4.0  APPLICATION  TO  THE  REENTRY  CONE 


The  techniques  just  outlined  can  be  used  in  the  problem  of  model 
simulation  for  the  reentry  cone.  Here,  one  can  take  advantage  of  the 
possibility  of  neglecting  various  terms  in  appropriate  regimes.  For 
our  purposes,  we  shall  restrict  our  application,  implicitly,  to 
simulation  in  a ballistic  range.  In  addition,  we  assume  (incorrectly, 
but  with  calculable  error)  that  there  is  exact  duplication  of  tempera- 
ture; i.e.,  Tm  = TTegt. 


4:1  GENERAL  CONSIDERATIONS 


If  flight  velocity  is  duplicated,  then  collision  efficiencies  and 

temperature  are  also  essentially  duplicated.  Furthermore,  at  high 

altitudes,  where  viscous  effects  are  minor  and  three-body  effects  are 

negligible  compared  with  two-body  effects,  simulation  is  accomplished 

(2) 

by  Mach  number  M and  by  D or  p H scaling.  Conversely,  at  low 

(3)  2 

altitudes  where  three-body  effects  might  dominate,  D or  p#  t scaling 
would  be  more  appropriate;  however,  the  additional  importance  of 
viscous  and  unsteady  influence  on  the  flow  field  serves  (as  has  been 
shown)  to  complicate  considerably  the  modeling.  Figure  2 qualita- 
tively suggests  (Ref.  5)  these  regimes  for  slender  cone  flow. 
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10"5  10"4  10"3  10"2  10" 1 1 
Ambient  Air  Density,  P^/Pq 


Figure  2.  Approximate  scaling  regimes  (from  Ref.  5). 

Other  altitude  related  complications  exist  with  regard  to  flow 
observables  which  add  to  the  modeling  difficulties.  Figure  3 shows, 
for  comparison,  a cone  and  a sphere  of  equal  drag  area  (hence  producing 
approximately  equal  numbers  of  electrons).  The  variety  of  behaviors 
for  the  slender  cone  is  much  greater  than  that  for  the  sphere.  For 
example,  there  are  four  regimes  in  so  far  as  wake  effects  are  concerned: 

(1)  high  altitude  - laminar  boundary  layer,  modest  electron  density; 

(2)  medium  altitude  - laminar  boundary  layer,  and  moderate  electron 
density;  (3)  medium  altitude  - laminar  boundary  layer,  and  high  elec- 
tron density;  (4)  medium  to  low  altitude  - turbulent  boundary  layer,  and 
high  electron  density.  Figure  4 qualitatively  illustrates  these  regimes 
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in  terms  of  wake  velocities  from  backscattering  radar  signals  and  shows 
a comparison  with  an  alternative  history  for  another  reentry  vehicle 
shape. 


Figure  3.  A schematic  for  shape  influences  on  wake  effects, 
(Cq  A)cone  — (Co  A)iphgra  x Mb  » 1. 


Using  Fig.  3,  it  is  possible  to  explore  more  explicitly  the  applica- 
tion of  binary  scaling  to  hypersonic  conical  flow.  For  example,  in  the 
inviscid  shock  layer,  if  local  flow-field  temperatures  can  be  held 
constant  by  preserving  M (sin  0),  then  not  only  is  useful  but  p^A/w 
is  a good  scaling  parameter  as  well  (Ref . 6) . This  allows  a direct 
scaling  of  shock  layer  radiation.  Notice,  however,  that  modeling  the 
radar  backscattering  still  presents  problems  since  the  plasma  frequen- 
cies depend  on  the  square  root  of  electron  density,  and  the  appearance 
of  critical  density  effects  would  not  be  modeled. 
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The  period  of  constant  velxity  would  correspond  to  underdense 
scattering  and  minimal  deceleration  of  the  reentry  vehicle. 


Altitude,  km ► Altitude,  km ► 

a.  A sharp-nosed  body  b.  A blunt-nosed  body 

Figure  4.  A schematic  for  shape  influences  on  altitude  histories; 

wake  velocity  measurements  at  a fixed  distance  from  the  nose. 

In  the  boundary  layer,  binary  scaling  should  still  dominate  at 
high  altitudes  if  there  are  no  unique  surface  effects.  This  requires, 
specifically,  a non- turbulent,  non-ablative  boundary  layer  so  that  Re 
and  St  can  be  ignored.  However,  a contaminating  boundary  layer  (in  order 
to  be  scaled)  would  require  the  scaling  of  viscous  and  heat  transfer 
phenomena.  Thus,  the  modeling  of  flow  with  binary  scaling  at  high 
altitudes  where  ablation  and/or  seeding  are  present  requires 
Re,  and  St  similarity. 

With  regard  to  the  reentry  wake,  there  seem  to  be  three  distinct 
regions;  (1)  base  flow,  between  the  cone  and  the  neck;  (2)  far-wake 
flow,  where  asymtoptic  laws  are  obeyed;  and  (3)  the  region  intermediate 
between  (1)  and  (2)  above.  Even  at  high  altitudes,  it  is  likely  that 
three-body  and  two-body  processes  are  at  least  competitive  in  regions 
(1)  and  (2).  However,  the  proper  scaling  for  the  electron  density 
decay  in  region  (3)  is  probably  characterized  by  Fig.  2.  Thus,  at 
high  altitudes,  the  rate  of  dissipation  in  radiation  and  (for  under- 
dense scattering)  in  radar  backscattering  should  follow  binary  scaling 
even  though  the  absolute  levels  of  electron  densities  may  not  scale. 
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There  are  several  additional  comments  which  are  appropriate. 

First,  the  usual  choice  for  a length  scale  would  be  £ ^ 'Jc^A.  This 

allows  for  the  possibility  of  investigating  modeling  effects  (assuming 
2 

CD  v sin  8)  in  the  scaling  of  wake  phenomena  which  could  be  attributable 
to  changes  in  the  cone's  angle.  However,  the  use  of  this  parameter 
implies  a truly  sharp  cone  (Rn/Rs  <<  1).  Secondly,  it  is  generally 
conceded  that  diffusion  effects,  molecular  and  turbulent,  are  not  prop- 
erly treated  in  the  similarity  analyses  (Ref.  7).  In  particular,  for 
example,  intense  small-scale  turbulence  could  produce  relaxation  effects 
which  would  not  appear  in  this  approach.  Furthermore,  the  flow  sequence 
at  the  base  of  the  cone  from  stagnant  flow  to  neck  flow  to  laminar  tran- 
sition to  turbulence  is  probably  inadequately  treated  in  this  approach; 
this  is  an  important  limitation  since  the  position  of  the  neck  plays  a 
crucial  role  in  determining  the  absolute  level  of  electron  density  in 
the  wake  and,  therefore,  the  quality  of  the  radar  backscattering. 


4.2  AN  ILLUSTRATION 

Now  let's  apply  these  results  to  the  problem  of  a prototype  high- 
altitude  cone  (50  km  < h < 90  km)  with  0 “ 10  deg  and  w = 6.1  km/sec.  It 
is  to  be  modeled  by  a ballistic  range  conical  projectile  with  8 = 10  deg 
but  length  scale  off  by  a factor  of  approximately  7.5.  Further,  allow 
for  a generator  inside  the  projectile  injecting  hot  gas  into  the 
boundary-layer  flow  at  the  base.  The  test  conditions  should  be  (with 

Too  “ TTest^  * 

(1)  Mach  number  duplication,  or  wTegt  ^ 6 km/sec 

£ 

(2)  Binary  scaling;  for  pJL  = PTest  then 

pTest  ' 7‘5  p- 
This  latter  condition  gives 


8.1  x 10  6 - p_  . - 8,1  x 10  8 g/cm3  or 
lest 


0.04  - PTegt  - 4 mm  Hg 
for  the  altitude  range  specified. 
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The  ballistic  range  observables  will  then  have  the  following  rela- 
tionship to  full  scale: 

1 . Radiation  from  the  inviscid  shock  layer  will  scale  with  -JC^A 
since  there  will  be  Mach  number  duplication. 

2.  In  the  absence  of  a firm  prediction  on  critical  electron 
density,  the  radar  backscattering  from  the  inviscid  shock 
layer  will  probably  not  scale. 

3.  In  the  absence  of  Re  and  St  similarity,  the  injected  hot  gas 
will  not  produce  a scalable  effect  on  absolute  electron  density 
or  radiation  since  binary  scaling  is  presumed  applicable  to  the 
prototype  problem. 

4.  Radiation  and  radar  backscattering  in  the  base  region  and  in 
the  far  wake  will  not  scale. 

5.  Relative  radiation  decay  will  scale  in  the  wake  for  interme- 
diate regions. 

6.  Relative  decay  in  radar  backscattering  will  scale,  assuming 
that  the  scattering  is  underdense. 

7.  Ionization  trail  lengths  will  not  scale  because  of  effect 

. . (3)  above,  viz,  the  nonscaling  of  absolute  electron  densities 
and  radiation. 

5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

, When  the  most  general  approach  is  taken,  the  cone  similarity 
problem  can  be  defined  in  a straightforward  fashion.  Some  limiting 
cases  are  especially  amenable  to  analysis  even  though  these  cases 
might  not  be  of  present  day  relevance.  The  possibility  of  usefully 
applying  the  techniques  of  similitude  to  carefully  selected  specific 
problems  has  been  shown  and  should  be  evaluated  since,  for  example, 
modeling  and  simulation  are  important  issues  in  threat-decoy  discrim- 
ination. It  is  recommended,  therefore,  that  experimental  tests  of 
predictions  of  the  sort  which  have  been  illustrated  be  conducted.  In 
addition,  these  procedures  should  be  extended  to  some  of  the  several 
blunt  body  types. 
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